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Abstract 

A quantum Markovian master equation is derived to describe the current 

noise in resonant tunnelling devices. This equation includes both incoherent 

and coherent quantum tunnelling processes. We show how to obtain the 

population master equation by adiabatic elimination of quantum coherences in 

the presence of elastic scattering. We calculate the noise spectrum for a double 

well device and predict sub-shot noise statistics for strong tunnelling between 

the wells. The method is an alternative to Greeen's functions methods, and 

population master equations for very small coherently coupled quantum dots. 
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I. INTRODUCTION 



Quantum features of conductance in mesoscopic electronics is currently a major theo- 
retical and experimental research interest in condensed matter physics ffl. Developments 
are driven by two complementary imperatives. Firstly a technological trend to fabricate 
devices on smaller and smaller scales is rapidly approaching the point where quantum ef- 
fects will become a problem unless explicit attempts to exploit quantum features are made. 
Quantum tunnelling can lead to undesired coupling between fabricated structures. On the 
other hand, tunnelling offers the possibility of very fast switching times. Secondly, the new 
devices require improvements to the theoretical description of electronic transport in a low 
temperature, high mobility regime. Small devices with very long coherence times can be 
dominated by coherent quantum effects. It is becoming increasingly clear that intrinsic 
quantum fluctuations play an important role at low temperatures || 

Current noise in resonant tunnelling devices (RTD) provides a path to understanding 
noise in the deep quantum domain. In a biased RTD one or more bound quantum states are 
coupled incoherently to two electron reservoirs maintained at different chemical potential. 
There are a number of experimental [§-|6| and theoretical |7Hl2f results. RTDs involve 
exchange of fermions between the reservoirs, and the bound states. We propose in this 
paper an approach to such devices based on quantum markov master equations [[□[]. Such 
an approach to quantum noise in nonequlibrium systems has been used with great success in 
quantum optics. This provides an alternative approach to the conventional Green's functions 
methods, and offers additional physical insights. For example it enables one to deal with 
coherent coupling between adjacent well states which couples off-diagonal elements of the 
density matrix in the occupation number basis and cannot be described by population 
master equations. Such coupling can occur in strongly coupled quantum dots, as in the 
recent experiments of Blick et al |14| and Oosterkamp et al |15| . 

If the strength of this coherent coupling dominates the time scales of elastic and inelastic 
relaxation, a population master equation cannot describe the system. Coherently coupled 
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Nanostructures are likely to become increasingly important and thus there is considerable 
motivation to develop theoretical schemes that go beyond population rate equations. 

In the first part of this paper we derive the operator master equation to describe a bound 
electronic system coupled incoherently to two reservoirs. We then apply this equation to 
calculate the current two-time correlation function for a single well, with a single bound 
state. This model of course can equally well be treated by a population master equation 
approach as in the approach of Carlos Egues et al. [16|, but we rederive the known results 
simply to display the method in a familiar context. In section III we apply our methods 
to treat the case of coherent coupling between the bound states of adjacent wells. In this 
case our approach yields results that go beyond the traditional population master equation 
approach. We derive the current spectrum in the device and demonstrate new features that 
arise precisely because of the coherent coupling between the two wells. To make contact 
with previous work we show that in the limit when elastic scattering dominates the coherent 
coupling, a population rate equation may be derived that is equivalent, in the appropriate 



limit, to that obtained by Carlos Egues et al. JTB 



II. THE MASTER EQUATION 

We begin with the derivation of the master equation for a single quantum tunnelling 
channel connecting two reservoirs under external bias. This system is quite adequately 
described by other methods, including population master equations. However we treat it 
here simply to demonstrate our approach in a familiar setting. Therefore our results are 
not new and could equally well be obtained by other methods. This is not the case for the 
coherently coupled double well system we discuss next. The hamiltonian describing this 
process is given by |J 

H = £ k a k a k + £ c c ] c + e p b l b P (!) 
k V 

+ J2( T Ek^a k + T* k a{c) 
k 
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+ J2(Tc P blc + T£ p c%) 
v 

where ak(a\.), c(c^) and b p (tf ) are the annihilation (creation) operators of electrons in the 
emitter (E) reservoir, in the central quantum well and in the collector (C) reservoir re- 
spectively. The energy of the bound state without bias is Eq which under bias becomes 
£ c = — ae V where a is a structure dependent coefficient. The single particle energies in 
the emitter and collector are respectively, = k 2 /2m and = p 2 /2m — eV. The energy 
reference is at the bottom of the conduction band of the emitter reservoir. 

The fourth and fifth terms in the Hamiltonian describe the coupling between the qua- 
sibound electrons in the well and the electrons in the reservoir. The tunnelling coefficients 
TEk-,Tcp depend on the barrier profile and the bias voltage. We will assume that at all 
times the two reservoirs remain in thermal equilibrium, with chemical potentials hc->He, 
with he — He = eV^, despite the tunnelling of electrons. This is one of the key defining 
characteristics of a reservoir. It assumes in effect that two very different time scales describe 
the dynamics of the reservoirs and the quaisbound quantum state in the well. 

In the interaction picture the Hamiltonian may be written as 

2 

Hj{t) = h^^T^e^ 01 + cT\(t)e- iuJot ) (2) 
where the bound state frequency is c^o — £c/^ and the reservoir operators are given by 

T 1 (t) = J2T Ek a k e-^ t (3) 

k 

r 2 (t) = E^ P e~^ (4) 

p 

where 

tf = £ (5) 

< = | (6) 

We now obtain an equation of motion for the density operator of the bound state, p(t), 
in the well following the standard method based on second order perturbation theory and 



tracing over reservoir states [17 . Thus we need 



^ = -jpfifaTrdHft), [H^lpn^pit)}} (7) 

where p R is the thermal equilibrium state of the two reservoirs, and Tr R denotes a trace 
over the reservoir variables. Note that the factorisation of the well state and reservoir states 
has been assumed. This is reasonable if the well state and the reservoir states are initially 
uncorrelated and provided there is a wide separation in the relaxation time scales of the well 
state and the reservoirs. The only non zero correlation functions we need to compute are 



iEiit) = / < ^i(r t 1 (t 1 )r 1 (t 1 ))e-^(^) (8) 

Jo 

i E 2(t) = / < rft 1 (r 1 (t 1 )r t 1 (t 1 ))e-^(^) (9) 

Jo 

ia(t) = / < rft 1 (rt(t 1 )r 2 (t 1 )) e ^^) (10) 

JO 

i C2 (t) = f : ^ 1 (r 2 (t 1 )r|(t 1 )) e -^(*-* 1 ) (ii) 

Jo 



In order to illustrate the important physical approximations required in deriving the master 
equation we will now explicitly evaluate the first of these correlation functions. 

Using the definition of the reservoir operators and the assumed thermal Fermi distribu- 
tion of the electrons in the emitter we find 

W*) = J2^Ek\T Ek \ 2 fdhe^-^ (12) 

As the reservoir is a large system by definition we can replace the sum over k by an integral 
to obtain 

fOO rji i FT , 

I E i(t)= ^p(luMlu)\T e (lu)\ 2 dre**-*** (13) 
Jo 2n Jo 

where we have changed the variable of time integration. The dominant term in the frequency 
integration will come from frequencies near ujq as the time integration is significant at that 
point. We assume that the bias is such that the quasibound state is well below the Fermi 
level in the emitter. This implies that near u = uo, the average occupation of the reservoir 
state is very close to unity. This is an effective low temperature approximation. Now we 
make the first Markov approximation. We assume that the function p(a>)n(a;)|T E (a;)| 2 is 
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slowly varying around u = Uo, and thus the frequency integration will lead to a function 
which is a rapidly decaying function of time compared to dynamical time scales for the 
quasibound state. This implies that on time scales of interest in an experiment we can 
extend the upper limit of the time integration to infinity as a good approximation. In that 
case I ei becomes time independent and may be approximated by 

roc A 11 

^p(uj)\T e \ 2 S(uj ) 
Jo 2 

= IlM (14) 

which defines the effective rate 7l of injection of electrons from the left reservoir (the emitter) 
into the quasibound state of the well. This rate will have a complicated dependence on the 
bias voltage through both u and the coupling coefficients \T e {oj)\- In this paper we do not 
address this issue. We simply seek the noise properties as a function of the rate constants. 

Evaluating all the other correlation functions under similar assumptions we find that the 
quantum master equation for the density operator representing the well-state in interaction 
picture is given by 

dp r 

Tt =cp 

= "t^(2c^ pc — cc 1 p — pcc^) 

+ -^(2cpc t -c f cp-pJc), (15) 

7l and 7_r are constants determining the rate of injection of electrons from emitter into the 
well and from the well into the collector respectively. The rate constants can be determined 
by a self-consistent band calculation involving the bias voltage. Two Poisson processes 
shown in the master equations: the injection of electrons into the well described by the first 
term in the right hand and the emission of electrons out of the well by the second term, are 
conditioned by the rates 7 i (cc l ~)(t) and 7#(c l ~c)(t): 

E (dN E (t)) = -y L (cc*)dt; E (dN c (t)) = lR (Jc)dt (16) 

where the average is taken with respect to the well state at any time t. The master equation 



T5| is diagonal in the occupation number representation. The mean occupation number 
n = Tr (ctcp(t)^ can therefore be determined easily from the rate equation 

= 7i(l - n) - -y R n (17) 

However the occupation number of the well states is not directly measured in current exper- 
iments. The current noise is a fluctuation in classical stochastic processes. It is measured in 
the relatively high temperature reservoirs of the leads, well away from the well state, and the 
strong, fast electron-electron interactions in the reservoir establish the classical level of the 
observed variable. It is however conditioned on the underlying quantum stochastic process 
in the well, which is described by the master equation. We thus have the familiar problem 
of connecting the observed classical stochastic process to the quantum source of information 
in an open quantum system. In this problem we proceed as follows. The current pulse in 
the emitter and collector may be determined from the Ramo-Shockley theorem fl8| . For a 
symmetric geometry this takes the form, 

t(t)dt=^(dN E (t)+dN c (t)) (18) 

The connection to the quantum source is then made by equations (|I(J). Using Eq(16), the 
average current is given by E{i{t)) = 7^(1 —n) —■ynn. In the steady state this is ioo = nL ^ R 
where 7 = 7l + 7_r and the subscript oo indicates the steady state. 

The fluctuations in the observed current, i(t) are quantified by the two-time correlation 
function: 

G(r) = ^ oo 5(r) + (I(t),I(t + r)y c i (19) 

To relate these classical averages to the fundamental quantum processes occurring in the 
well we apply the theory of open quantum system |0J to the present system and calculate 



the following correlation components with r > [20]: 



E (dN E (t + T)dN E (t)) = -i\Tr (cc 1 e Cl c 1 Poo c) dt 2 (20) 
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E (dN c (t + T)dN c (t)) = 1r Tt (c ] ce CT c Poo c ] ) dt 2 (21) 
E {dN E {t + r)dN c {t)) = lLlR Tr (cc 1 e Cl c Poo J) dt 2 (22) 

E (dN c (t + r)dN E (t)) = lRlL Tr (Vce^cW) dt 2 (23) 

Calculating the above correlation components using master equation with corresponding ini- 
tial conditions, and substituting them together with the shot-noise component into equation 
U yield: 



G{r) = — 5{t) + — I 1 —)l e7n (24) 



Thus the spectral density of current fluctuation in frequency domain is given by 



S(u) = 2^°° G(t) (e iUT + e~ iujT ) dr 



The current Fano factor -F(^) is defined as the ratio of current noise density over the full 
shot noise density, and for low frequencies (7 ^> u): 



The shot noise is suppressed and reaches the minimum of 50 % in a symmetric structure with 
7l = 7.R- The result is same as those derived by Chen and Ting [[| using non-equilibrium 
Green's function method. The result can also be obtained by a classical master equation 
calculation [JCTT]. However the classical master equation cannot be used to treat the case of 
coherent coupling in a double well system discussed below. 

The suppression of fluctuations at low frequency is due to the exclusion principle in the 
well state, reflected in the master equation by the appearance of the anti-commuting field 
operators. No electron can tunnel onto the well if an electron is already there. We need to 
wait a time of the order of 7 _1 for the electron to tunnel back out into the collector. Thus 
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strong anti-correlations are established in the two fundamental Poisson processes, dNi(t). If 
the tunnelling particles were bosons, the well could accumulate a large number of particles, 
enhancing the probability for emission into the collector. This would lead to a rapid bunching 
of emission events into the collector and a super-shot noise current would result. At high 
frequencies, we are looking at fast processes in which an electron tunnels into the well and 
immediately tunnels out. The Fano factor at high frequencies is 0.5 due to the assumed 
form of Ramo-Shockley theorem. 

III. NOISE PROPERTIES OF A COHERENTLY COUPLED DOUBLE WELL 

STRUCTURE 

We now apply our approach to a triple barrier and double quantum well involving elastic 
scattering within the wells and coherent coupling between the wells. The main procedures 
are parallel to those in single well case but now involve off-diagonal elements of the density 
matrix. The master equation is, 



where Ci(c{) and c 2 (d>) are annihilation (creation) operator of electron in the left and right 
quantum well respectively, rji is the rate of elastic scattering (electron-phonon for example) 
in the ith well, and Q is the coherent coupling rate between the two well states. The 
irreversible term describing the elastic scattering is derived in much the same way as the 
inelastic tunnelling terms that describe electrons entering and leaving the device, with one 
additional assumption. To get a Markov master equation for number conserving scattering 
events we must assume that the temperature of the bath describing such processes is high 





c[d,p 



-V2 |c 2 c 2 , [C 2 C2, P\ 

-in (4c 2 + 4ci),p 



(27) 
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enough that the bath states are well away from the ground states. This is not a very 
restrictive assumption for realistic devices at milliKelvin temperatures. The deviations that 
can result for very low temperatures are described in Gardiner [HJ The derivation of the 
scattering term in the master equation (27) is detailed in the appendix. 

The last term in this equation represents a coherent coupling between the two wells 
and causes a single electron to periodically tunnel backward and forward between the two 
wells, until it is eventually lost through the final barrier. Recently Blick et al fl4| , have made 



measurements on a structure that can be roughly approximated by our model. As they point 
out this device exhibits a new feature, in that a single electron can be in a superposition 
state between the two wells and is thus like an artificial molecule. We first derive the noise 
features in the presence of this coherent coupling. We will then show that, in the limit of 
strong elastic scattering, rji >> Q , the system can be described in terms of population 
rate equations that have been extensively used in the past. 
The steady state current is easily found to be given by 

~ 7 2 + 2 leVe + AW (28) 

for a symmetric system, 7^ = 7^ = 7e ; r\\ = r] 2 = r\ e . The appropriate correlation functions 
may be evaluated to give, 



A+A 
1 

4A 



2 



<I(f),I(t + r)) M =(^fJ {4fi 2 

/ + e A + T + /_e A - T ]} ( 29 ) 



where 



f± = (7e - Ve ± A) (A ± tj e ) ( 7e + tj 6 + A) (30) 

and A = \Jrjl — 4f2 2 ; A± = — 7e — r\ e ± A. The noise spectra are derived in two cases. In 
easel: r/ 2 > 4f2 2 , when the elastic scattering rate, rji, is higher than the coherent coupling 
rate between the well states, the current noise spectrum is 
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(r) e + A)( 7e -77 e + A) 



, (-T] e + A)(7 e — TJe — A) „ , ■ 

+ (-^-r^-A) 2 ^ 2 Jl (31) 
The current Fano factor against normalised frequency is plotted in Fig. 1 where the 
spectrum shows Lorentzian feature. In case2: if e < 4Q 2 , the opposite situation, when 
coherent coupling is much stronger than elastic scattering, the noise spectrum 



( 7e -77 e + 2A)(77 e + iA) 

(~7e - Ve + ^A) 2 + CJ 2 



(32) 



where A = y4fi 2 — r^, is symmetrically double peaked about the free particle frequency 

{uj = 0) as shown in the Fig. 2. A comparison of these two quantum processes is shown in 

Fig. 3. When elastic scattering increases, the Fano factor increases. Increasing the coherent 

coupling results in noise suppression and the double peak feature are more significant as the 

two well coupling and the quantum correlations are stronger. Further when elastic scattering 

is extremely weak: r) e — > 0, and coherent coupling is strong: Q ^> j e , the steady state current 

ioo — > — approaches the single well case as expected in this limit. A significant outcome is 
2 

that the best noise reduction at low frequency when r/ e = reaches 0.22. 

The coherent tunnelling between the two wells has a strong effect on the noise char- 
acteristics. Electrons are periodically transferred between the two wells at the tunnelling 
frequency. If an electron from the emitter is injected into the first well, no further electrons 
can enter this well until this electron is removed, which takes place on a time scale deter- 
mined by Q^ 1 . Thus at frequencies smaller than Q, noise is suppressed by the exclusion 
principle, just as for the single well case. At the tunnel frequency however we expect the 
noise to increase, as electrons injected into the first well are quickly cycled to the second 
well, where they can incoherently escape to the collector. This explains the two peaked 
structure of the noise power spectrum. In the case of large Q however, coherent coupling 
dominates. In that case if an electron tunnels into the first well it periodically returns to 
that well at a frequency of 2Q. To see this it is sufficient to note that the two levels which 
are degenerate in the absence of tunnelling become split into symmetric and anti symmetric 
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combinations, separated in energy by 2M1 A state initially localised in one well can then 
be written as a linear combination of the two new eigenstates. The phase difference in the 
superposition rotates through n at the frequency Q which leads to a state localised in the 
other well. This is just the standard description of tunnelling in a two state system. The 
periodic of return of the electron to the first well suppresses another electron from entering 
the well. Thus at large values of fl we expect noise suppression to occur at uj = 2Q. This 
behaviour is indeed seen in figure4. 

We now show that in the limit of strong elastic scattering rji » Q (Case 1 above), a 
population master equation can be derived that describes a classical sequential tunnelling 
structure. The sequential model is traditionally formulated in terms of a classical master 
equation for the occupation probabilities of each well. In our case, we have restricted the 
discussion to a single bound state in each well and thus the maximum population in each 
well is unity. However we can derive an equivalent classical master equation to describe 
sequential tunnelling even in this case. 

Our method is an extension of adiabatic methods used in quantum optics to obtain rate 
equations. We assume that the off-diagonal elements of the double well density operator are 
rapidly damped due to the elastic scattering rates rji. The off-diagonal elements are then 
assumed to relax almost instantaneously to their steady state values and adiabatically follow 
the more slowly changing diagonal matrix elements. 

From Eq(27) we find the following equations of motion for the matrix elements in the 
occupation number basis for each well, 



d 

< n 1 n2\p\mim2 > 



dt 



= < 



- y[( n i + Wmfl + ( m i + l)<*mi,o] - y (™2<W + m 2 5 m2jl ) 



-r)i( n i^n u i - 2niTOi5 nij i5 mi) i + m\8 mul ) 

-772(^^2,1 - 2n 2 m 2 <5 n2i i<5 m2i i + m\b m ^ x ) 
+lL6 ni ,i5 mi ,iy/n 1 m 1 < n x - 1, n 2 \p\m 1 - 1, m 2 > 



> < nin 2 |p|mim 2 > 



+7^n 2 ,o ( W(-l) ni+mi + l)(m 2 + 1) < ni,n 2 + l\p\m 1 ,m 2 + 1 > 
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-iQ{5 nul 5 n2fi (-l) ni - 1 ^n 1 (n 2 + 1) < n x - l,n 2 + l\p\m 1 m 2 > 

+5 ni ,oSn 2 , i(~^) ni \j i n i + l ) n 2 < ™i + 1,^2 - l\p\m 1 m 2 > 
S mifi 5 m2A ( y -l) mi ^( y m 1 + l)m 2 < n 1 ,n 2 \p\m 1 + l,m 2 - 1 > 
-5 mi ,i^m2,o(-l) mi_1 \/(^2 + l)mi < n 1 ,n 2 \p\m 1 - l,m 2 + 1 >} 



where rii, n 2 refer to the occupation number of the first and second wells respectively. 
Note that the diagonal matrix elements represent the occupation probabilities of each 
well, 

P{n u n 2 ,t) = (ni,n 2 |p(£)K,n 2 ) ( 33 ) 

The diagonal matrix elements then obey the equation, 
d 

— < nxn 2 \p\nxn 2 > 
at 

= [-7l(™i + l)^ni,o - lRn 2 bn 2 ,i} < n x n 2 \p\n x n 2 > 
+S ni ,ilLni < ri! - l,n 2 \p\n 1 -l,n 2 > 

+Sn 2 ,0lR(n 2 + 1) < 711,712 + l|p|ni,7l 2 + 1 > 

+i(-l) ni ^{S nijl 5 n2fi ^n 1 ( y n 2 + 1)[< n x -l,n 2 + l|p|nin 2 > - < n x , n 2 \p\n x - 1, n 2 + 1 > 
+<5m,o<WyVi + 1W< ni,n 2 |p|ni + l,n 2 - 1 > - < n x + l,n 2 - l|p|nin 2 >]} 
we now define the off-diagonal matrix elements as 



Y 1 =< ni,n 2 \p\rii + l,n 2 - 1 > 

Y 2 =< m, n 2 \p\ni - 1, n 2 + 1 > 

Y 3 =< ni - l,n 2 + l\p\ni + l,n 2 — 1 > 

Therefore, the population equation we are interested is, 

-^P(n!,n 2 ,t) 

= [-7l(«i + !)^ni,o - 7i?^25„ 2 ,i)]P(ni, n 2 , t) 
+<5ni,l7L"'l-P(™l ~~ n 2 ? + Sn 2 fllR{n 2 + l)P(ni, n 2 + 1, t) 
-2fi(-l) ni [<5„ li0 5„ 2 ,i v /(^i + l)n 2 \vciY 1 - 6 nul 5 n2fi ^/ j n x {n 2 + 1) Im7 2 ] 

13 



Note that the elastic scattering rates, 771,772, do not directly enter this equation. This is 
because elastic scattering does not change the occupation of the well states but does disrupt 
the phase coherence between the wave functions in the wells. This will lead to a decay of 
the relevant off-diagonal matrix elements, which obey the equations, 



d , « d 

H *i(0 = H < n i n 2\P\ni + l,n 2 - 1 > 
at dt 

= { 2"("1 + !)^ni,0 - — ^25n 2 ,l - 77l[7Z^ ni; i + (m + l) 2 5„ 1+ i,i)] - r] 2 n 2 2 5 n ^ l}*l(*) 



^"l-i 



tt5 nit i5 n2t o^ni(n 2 + l)Y 3 (t) 



-i(-l) n ^5 nifi 5 n2 ,i V(m + lKP(ni + l,n 2 - 1,0 



+?(-l)" 1 l¥ ni+ i i i5 n2 _i ii0 y / (ni + l)n 2 P(ni, n 2 , t) 
d d 

-j: Y 2{t) = -77 < 7li7l 2 p Tli - l,7l 2 + 1 > 

= {- + 1 ) <J m,o + wi<J ni _i,o] - ^[n 2 5n 2 ,i + (n 2 + 1)^2+1,1] 
-VinjSn^i - rj 2 [nl5 n2:1 + (n 2 + l) 2 <5 ri2+ i 5 i]}y 2 (^) 

-i(-l) ni n5„ 1 ,o5n 2 ,iV( ni + lK y 3*(0 
-^-l)" 1 - 1 ^,!^^^^ + l)P(ni - l,n 2 + 1,0 
+i(-l) ni - 1 tt5 ni - lfi 5 ri2+1:tl ^n 1 (n 2 + l)P(n u n 2 , t) 

d d 

37^3(0 = — < Tli - 1,tz 2 + l|p|ni + l,n 2 - 1 > 
dt dt 

= [ — - yh + l)<5na+i,i - + l) 2 ^m+i,i - ^2(^-2 + l) 2 S n2 +i,i]Y 3 (t) 



-i{-l) n ^5 n ^ lfi 5 n2+1 ^ ni (n 2 + 1)^(0 + i(-l) ni ^„ 1+ i,i5„ 2 -i„oV(«i + 1)^(0- 

To proceed we solve the equations for Y 1 ,Y 2 ,Y 3 in the steady state, assuming that the 
diagonal matrix elements are time constant in time over the lifetime of the of-diagonal ma- 
trix elements. This is the adiabatic approximation. These steady state values are then 
substituted back into the equation for the diagonal matrix elements to obtain a classical 
jump process master equation to describe sequential tunnelling. The algebra is tedious, so 
we will not give details. The result is 
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d ^ . 

—P{n 1 ,n 2 ,t) 
at 

+*ni,iT£«i- p («i _ ra 2, *) + 5n 2 ,o7i?(^2 + l)P{n 1 ,n 2 + 1, t) 
-2fi 2 {5 ni , 5 n2 ,i(n 1 + l)n 2 (Q22Q33 ~ a23Q32) [P(m + l,n 2 - 1,0) - P(n 1; n 2 ,0)] 
+5 ni)1 5 n2|Q n 1 (n 2 + ^ (^33-^^31) ^ _ 1;n2 + 1>0) _ P(ni)n2)0 )]}. 



where D is given by 

D = a n (a 22 a 33 - 0236132) - a 13 a 22 a 3 i (34) 

with 

a n = — — (ni + l)5 nifi - —n 2 8 n2tl - r]i[n\8 nu i + (rai + l) 2 5 nil o) - rj 2 n 2 2 5 n ^i 



a i3 = (-l) ni fi(y ni ,l*na,0V n l( n 2 + 1) = -%L 

«22 = - y[( n i + i^m.o + ni<y nij i] - ^p[n 2 <J„3,i + (n 2 + 1)5„ 2 , ] 



-^l^f^ni.i - r} 2 [nl5 n2A + (n 2 + 1) 2 5 

"2,0j 



«23 = -(-l) ni ^ni,0^n2,iv( n l + !) n 2 = -«32 
«33 = 2 UlSni ' 1 ~ IT^ 2 + 1 )*"a. _ + 1 ) 25 ni,0 ~ ^2(^2 + l) 2 ^n 2 ,0 



In addition to the incoherent tunnelling of electrons between the wells and the external 
reservoirs, we now have incoherent (sequential) tunnelling between the two wells at rates 
determined by — . The form of this equation corresponds to the sequential tunnelling master 
equation obtained by Carlos Egues et al [jnj. We have thus shown that, in the limit of strong 
decoherence induced by elastic scattering of the bound states, a population master equation 
may describe sequential tunnelling in the device. This will be the appropriate limit in the 
case that r)i » Q. However future quantum nanostructures devices are likely to operate 
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in the opposite limit. In that case our method is ideally suited for determining the device 
characteristics. 

IV. SUMMARY 

We have shown how the quantum theory of open systems, formulated as a quantum 
stochastic process, enables the current noise spectrum to be calculated for mesoscopic tun- 
nelling devices. Our approach explicitly treats quantum noise properties of the charge carri- 
ers, and gives a simple intuitive picture to understand the results. As fabrication technology 
develops, quantum noise limited networks of coherent tunnelling devices, such as quantum 
dots and quantum point contacts, will become increasingly important. Such coherently 
coupled devices are essential for the implementation of a quantum computer, which must 
operate reversibly ||21|| . The full operator master equation methods we have demonstrated 
here provide a powerful description, including both diagonal and off-diagonal matrix ele- 
ments in the same equation. 

Our model does not treat the transverse unbound modes in the well of a realistic resonant 
tunnelling device. These can easily be incorporated by additional states in the well and 
additional jump process channels in the master equation. We have not done that here as 
we sought to derive the irreducible level of current noise in tunnelling devices. Our model 
may in fact apply to very tightly confined quantum dot structures which could conceivably 
be fabricated with a single bound well states at donor impurities. Further extensions of the 
model are also needed to treat the case where the well state is just below the Fermi level in 
the collector in which case the current noise acquire an additional temperature dependent 
classical component. These more general cases will be treated in a larger publication. 
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In this appendix we derive the master equation describing elastic scattering of the quasi- 
bound states of the well which cause a dephasing of the electron quasi-bound states but 
do not change their populations. The Hamiltonian for double-well system in Schroedinger 
picture is 

H = H + H T + H scat (Al) 

2 

H =J2 £ nClc n + £ k a W + Y £ p b l b P 
n=l k p 

+ u q a\a q + &{c[c 2 + eld) (A2) 
q 

H T = £(lWa fc + T* Ek a\c) + £(7b^c + T* p c%) (A3) 

k P 
2 

H scat = J2 C i°n T n (A4) 



71=1 



where 

r\ = YM q (al + a q ) (A5) 
q 

T 2 = YM q ((3l + (3 q ) (A6) 
q 

where a q , /3 q are Bose destruction operators describing independent reservoir oscillators. 
Note that each bound state in the well is coupled to an independent reservoir. This assumes 
that there are no correlations between well states due to the dephasing that takes place 
through elastic collisions. 

We will only consider here the derivation of the master equation arising from the elastic 
scattering of bound states and the harmonic oscillator reservoirs. The relevant part of the 



master equation is |13 



dp{t) 1 '•' 
dt h 2 Jo 



f drTr B [H scat {t), [H scat (r), p(r) ® p B \\ (A7) 
jo 



where p_B is the equilibrium state of the bath, and where Tr B means to trace over the bath 
variables. This equation is based on a second order expansion in the interaction energy 
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between the reservoir states and the bound states of the well. We have also assumed that 
the system and bath states are decorrelated very rapidly on the time scale of interest in the 
system, so that the bath remains close to its equilibrium state. The bath Hamiltonian is 

H B = ^ VqO>l a q ( A8 ) 

i 

It will only be necessary to consider one of the bath-well state coupling terms in the 
scattering Hamiltonian. The relevant part of the master equation in the interaction picture 
is 

q 

[c{c n J2 M q {a\e^ + a q e~^ T ), p(r) ® p B }}} (A9) 



We now define, 



E ee cic^M^aU^ + oye-%") (A10) 



F = cic n ]T M q {a\e^-^ + a^^) (All) 



Therefore 



= ~T2 f drTr B {EFp{r) ® p B - Ep(r) ® p B F 
at h Jo 

-Fp(r)® p B E + P (t)® p B FE} (A12) 

The state of the reservoirs is taken to be a thermal state at temperature T, thus 

Tr B {a q a q >p B ] = Tr B [a\a\p B ] = (A13) 
Tr B [a\a q ,p B \ = 5 qq> _ - (A14) 

The first term on the right hand of the equation (A12) 



o 



drTr B EFp(t)® p B 



11 



^i[(ul q ,~Ul q )t+Ul q T} 



18 



+(1 



1 



e E q /k B T _ I 



^-^ ]t -^ T] ](cic n ) 2 P 



E|M/[ sin(M (1 



+ 



fan +^n)(4c n ) 2 p 



e E q /k B T _ I 



UJ„ 



(A15) 



Where 



E|Mg|2 sin(M (1 



+ 



e E q /k B T _ l 



) 



(A16) 



2 cos(u q t) — 1 



(A17) 



Similarly, 



- f dTTr B Ep(r) ® p B F 
Jo 



fan ^njc^CnpC^Cn 



(A18) 



/"* dTTr B Fp{T) ® p B £ 

JO 



fan + ^„)4c rl p4c 



(A19) 



/* drTr B p{r) <g> p B FE 
Jo 

fan ~ i£n)p{cic n ) 2 



(A20) 



The coefficients r/ n , £„ appear to be time dependant, but under reasonable physical as- 
sumptions are time independant fl3fl . These assumptions are, firstly that t is assumed to 
be a time scale over which the system operators vary significantly. On this time scale bath 
correlation functions decay rapidly. Secondly, that the bath is at finite temperature and 
there is significant excitation above the reservoir ground state. Finally that the coupling 
constants M q are independent of q up to some large cut-off wave number. Under these 
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assumptions these coefficients can be evaluated in the limit of t — > oo. We refer the reader 



to reference [DJ for more details. 

The total contribution from the scattering term to the master equation is therefore by 
substituting equations (A15) and (A18 -A20) and coresponding termes for the second well 
into equation (A12 ): 

= J2Vn[cic n , [c{c n ,p]} + p X^n^P] ( A21 ) 
al n=l 11 n=l 

The effect of the term i£ is to add a small perturbation to the energy of each quasi-bound 
well state and is equivalent to the Lamb shift term in atomic physics. 
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FIGURES 

FIG. 1. The current Fano factor |^ versus normalised frequency to/j e in case 1. All param- 
eters are normalised by 7 e . The corresponding parameters (r] e /^ e , uj/^ e ) for the curves are from 
the top: Dotted:(0.5, 0.1); Dot-dash: (0.5, 0.2); Solid:(0.5, 0.24). 

FIG. 2. The current Fano factor versus normalised frequency (u}/j e in case 2. All parameters 
are normalised by j e . The corresponding parameters ij e /j e , fi/7 e ) for the curves are from the top: 
Dotted:(0, 0.2); Dot-dash: (0.4, 0.4); Solid:(0.5, 0.5). 

FIG. 3. The comparison of influences of the elastic scattering and the coherent tunnelling. The 
corresponding parameters (r] e /j e , O/Te ) f° r the curves are from the top: (a):(0.6, 0.1); (b):(0.0, 
0.2); (c):(0.2, 0.5); (d):(0, 0.645). 

FIG. 4. The current Fano-f actors versus normalized frequency for double well structures. The 
normalized parameter £l/j e for the curves are from top: 0.2 dash line , 0.645 dot-dash line, 5.0 
solid line. 
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